Mott insulating phases and magnetism of fermions in a double-well optical lattice 
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We theoretically investigate, using non-perturbative strong correlation techniques, Mott insulating 
phases and magnetic ordering of two-component fermions in a two-dimensional double-well optical 
lattice. At filling of two fermions per site, there are two types of Mott insulators, one of which 
is characterized by spin-1 antiferromagnetism below the Neel temperature. The super-exchange 
interaction in this system is induced by the interplay between the inter-band interaction and the 
spin degree of freedom. A great advantage of the double-well optical lattice is that the magnetic 
quantum phase diagram and the Neel temperature can be easily controlled by tuning the orbital 
energy splitting of the two-level system. Particularly, the Neel temperature can be one order of 
magnitude larger than that in standard optical lattices, facilitating the experimental search for 
magnetic ordering in optical lattice systems. 
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Introduction.- There are currently worldwide efforts in 
studying collective properties of cold atoms either in a 
single trap or in an optical lattice pQ. A central goal 
of these studies is to explore novel many-body quantum 
phases in both bosonic and fermionic systems. While 
both bosonic and fermionic Mott insulators have been 
realized in labor atories[2]-[5], the experimental search for 
magnetism in optical lattices is currently on-going. Most 
of these studies have been focusing on the single-band 
physics. For example, it is known that two-component 
fermions in the lowest band can be used to study spin-1/2 
antiferromagnetism[S] . 

A question naturally arises: Is it possible to realize 
multi-band magnetic systems using cold atoms in opti- 
cal lattices? Theoretical studies suggest, for example, 
exploring excited bands in optical lattices for searching 
novel magnetising |5], partly because of the enhanced 
tunnel coupling in excited bands [TO]. While there 
are currently experimental efforts along this direction to 
populate atoms in excited bands [TTHT2"] . whether one can 
overcome the finite-life time problem of atoms in excited 
bands still remains unclear. 

On the other hand, there is a crucial practical issue on 
the energy and time scales of atoms in optical lattices. 
Ordinary optical lattices are characterized by extremely 
small energy scales, which also lead to slow relaxation of 
lattice systems to equilibrium. For example, the tunnel- 
ing of the lowest band is about a few nano-Kelvin which 
corresponds to a time scale of a few tens of milliseconds. 
The energy scale associated with the super-exchange in- 
teraction t 2 /U is even smaller since U ^> t typically. As 
a result, the Neel temperature of antiferromagnetism in 
ordinary optical lattices is far too low for experimental 
observation. Meanwhile, it is also challenging for the sys- 
tem to reach equilibrium because of the long relaxation 
time. A scheme to enhance the relevant energy scales is 
therefore very desirable, particularly in the context of the 
experimental study of many-body magnetism in optical 



lattice systems. 

In this Rapid Communication, we theoretically study 
quantum magnetism of fermions in a double-well optical 
lattice. Instead of the usual spin-1/2 magnetic ordering 
in an ordinary optical lattice, the double-well effectively 
produces a spin-1 system. The associated magnetism is 
induced by the inter-band interaction between the lowest 
two bands, and is a ground state property. Moreover, the 
characteristic energy scale for observing magnetism can 
be enhanced by one order of magnitude compared with 
the spin-1/2 magnetism in ordinary optical lattices. As a 
result, the magnetism may be much easier to achieve and 
observe experimentally in the double- well optical lattice. 

A double-well optical lattice contains two potential 
wells, which are separated by a barrier, on each lattice 
site. Its unique advantage is that the band gap between 
the lowest two bands is tunable p~3 \ I14 j . When these two 
bands are very close to each other, interesting quantum 
many-body phenomena, which are completely absent in 
ordinary optical lattices, emerge [HI [IS]. We will see that 
the interplay between the orbital degree of freedom and 
the fermionic spin lies at the heart of the new physics 
reported here. Theoretically it is however challenging 
to study fermions with spin degrees of freedom in the 
presence of multiple bands. We employ the dynamical 
mean- field theory (DMFT) method [l^j, and study both 
the Mott insulating phases and magnetic properties of 
fermions in a double-well lattice. We show that at filling 
of two fermions per site, the Mott insulator developed 
in the system can be either the triplet (n s ,n p ) — (1, 1) 
states or the admixture u(2,0) — v(0, 2). For the for- 
mer case, antiferromagnetic order emerges in the spin-1 
channel, which should be experimentally observable. 

Model- We consider the Hamiltonian containing a 
tight-binding-band part and an on-site interaction part, 
characterizing the two lowest bands (labeled by s and 
p respectively) in a symmetric double-well lattice, H = 
-ffband + -ffint- In the real space, the band part can be 
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FIG. 1: (Color online) (a) Profile of the double-well potential 
(the green line) along the ^-direction. The s (red lines) and 
p (blue lines) orbitals are schematically shown. The hopping 
integrals and energies for s(p) orbitals, t sx (t px ) and e„(e p ), are 
indicated. Note that the potential is not double-welled along 
the y-direction (not shown), where the hopping integrals are 
t sy and t py correspondingly, (b) (n s ,n p ) = (1,1) eigenstates 
in atomic limit, including triplet states with energy Et and 
a singlet state with energy Es- (c) Linear combinations of 
(2, 0) and (0, 2) states as eigenstates for the Hamiltonian in 
atomic limit, whose energies are denoted by E±. See the text 
for details. 



written as 



-ffband = ^2 \(£ a - /U)4yr s <r,r + (s p - [l)p\ 
rcr 

~t~ ^ ^ ^ ^sx^^r &<T,v-\-x ^syS^a.r^&iT+y (1) 
rcr 

~T" tpxP^j.rPcr.r+x ^pyPcr : rPcr : r+y ~\~ h.C 



where s\ r (p\ r ) creates a fermion with spin a on the s(p) 
orbital of site r, e s and e p are the energies for s and p 
orbitals, and fi is the chemical potential. The hopping 
amplitude for s and p orbitals may differ in x and y direc- 
tions, thus we label them by t sx , t sy , t px , t py respectively. 
The interacting part of the Hamiltonian can be written 

as 



— ^ ' U s n s -\ t rn a \, t r + U p n p ^^ r n p ^^ 

V 

+ U sp (n 

~ U »P { S l,rP\,rPl,r S t,r + P\, r Pl,r S t,r S l,r + h - c -) 



(2) 



where n aiTjr = aj^ r oti a>r (a = s,p) is the number oper- 

fd 3 xW*(x) 



ator for orbital a at site r, U a 
denotes the intra-band interaction, while U s . 
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Jd 3 xW / s 2 (x)W p 2 (x) denotes the inter-band inter- 
action, where a s is the scattering length, M is the mass 
of the fermion, and W a (x) is the Wannier wave function 
for each band. The inter-orbital terms in Eq.([2]) char- 
acterized by U sp are referred as density-density, spin- 
exchange and pair-hopping interaction. This model is 
essentially the rotationally invariant Slater-Kanamori in- 
teraction widely studied in transitional metal oxides |17j . 
The main difference here is that the spin-exchange and 
pair-hopping are as strong as the inter-orbital density- 
density interaction. An important parameter which con- 
trols the multi-band physics is the energy level splitting 
between the two levels, defined as A = e p — e s . When 
A is small or intermediate, interactions between the two 
orbitals give rise to interesting phenomena, as we dis- 
cuss below. When A becomes very large the physics 
reduces to that of the single-band model. The Hamil- 
tonian in Eq. ^ has been previously considered for 
fcrmions at resonance in an ordinary optical lattice in 
one-dimension[18 . In our case, the reduced band gap 
makes the realization of a two-band system more practi- 
cal in current experiments. Moreover, higher dimension- 
ality of our system gives distinct physical phenomena not 
accessible in one dimension. 

We start from the atomic limit, where the tunneling 
terms are absent. We are interested in the states at fill- 
ing of two fermions per site, the schematics of which are 
shown in Fig. [ljb) and 1 (c) . When there is one fermion 
in each orbital, they form triplets which is denoted as 

(n s ,n p ) = (1,1), namely, pjsj |0), -j^ fpjsj + pjsj) I )' 

and |0) , with degenerate energy Et — 2(e s — fx) + A. 

The singlet state ^ (p\ s \ ~ -Pt s t) 1^ ^ as a hi§ ner en ~ 
ergy E s = 2{e s - fi) + A + 2U sp . E s > E T sim- 
ply because two fermions interact with each other by 
s-wave short-range interaction. In the spirit of the Hub- 
bard model, atoms with different spins repel with each 
other and atoms with the same spin do not interact. 
On the other hand, the two fermions can also form ad- 
mixtures u(2, 0) ± v(0, 2), as shown in Fig. |T](c), due 
to the pair-hopping interaction. The eigenenergies are 

±J(A+^Y^Ul. By 



E± = 2(e s -^)+A+ 1 ^ 



controlling A in the double-well lattice, E_ can be made 
either smaller or larger than Et- Throughout the paper 
we fix parameters U s — Ylt, U p — lit, U. 



sp 



12t, and 

vary A and the temperature T . Straightforward algebra 
reveals that at the critical value A c = At, E_ = Et- For 
latter use we note that the hopping integrals are chosen 

ctS ~^ sOu ' — " ^ s y ' — ~ '^fj'y ' — ~ ^ ^ ' — " 2^ • T^llO IOjT^G i St GUIS 

from the fact that p bands are spatially more extended 
along the a>direction. However, our solution to the lat- 
tice model as well as the physics therein does not depend 
on this particular set of parameters in any important way. 

When the hopping terms are switched on, we employ 
the single-site DMFT [TB] to solve the strongly-correlated 
interacting lattice fermion problem. The key approxima- 
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FIG. 2: (Color online) (a) Total occupancy (n s + n v ) versus 
chemical potential fi, calculated for three different values of 
A. The calculation is done without magnetic order at tem- 
perature T = 0.2t. e s = in this plot. The points (the circle, 
the square the diamond and the triangle) indicate the loca- 
tion at approximately the center of the gap for the line with 
corresponding color, where we study magnetic ordering, (b) 
The difference in occupancy n a — n p plotted at the same fi 
scale. 



tion is the neglecting of momentum dependence of the 
self-energy: £(fc,cj) — > £(w), which is solved iteratively 
from an auxiliary quantum impurity problem plus a self- 
consistency condition. We use the matrix representation 
of the continuous-time hybridization-expansion quantum 
Monte Carlo impurity solver [TH] prescribed specifically 
for multi-band interactions. This is a state-of-the-art 
highly demanding numerical solution of the strongly in- 
teracting multi-band lattice Hubbard model in the con- 
text of our double-well optical lattice system. 

Mott physics.- There are multiple choices to fill a single 
lattice site with two fermions, forming different types of 
Mott insulator. To distinguish them, we have calculated 
both n s + n p and ft g clS functions of // for different 

values of A, as shown in Fig. [2ja). A Mott insulating 
gap at filling two is evident for all cases, and there is no 
qualitative difference in the value of n s + n p between dif- 
ference cases. However, the difference in occupancy for 
the two orbitals (n s — n p ) shows distinct behaviors. At 
very large level splitting A = 4.5i, n s ^> n p . This is con- 
sistent with the analysis in the previous section for the 
atomic limit, where each lattice site is filled by the state 
it(2, 0) — v(0, 2) and u»». In contrast, at small A, e.g. 
as seen in Fig. [2jb) for A = 2.5i and A = 0, n s ss n p in 
the Mott insulating regime. This indicates that on each 
site the triplet states dominate the ground state. This is 
also consistent with the atomic limit where the energy of 
(1,1) states Et continuously decreases and eventually the 
(1,1) triplet becomes the ground state with decreasing A, 
as discussed in the previous section. The transition be- 



tween the two types of insulator is a crossover. Since in 
this paper we focus on properties at non-zero tempera- 
ture, we shall not discuss the nature of this transition at 
zero temperature. 

For A — > oo, the magnetism is manifestly absent and 
the ground state continuously connects to the trivial 
band insulator in the lowest band of an ordinary opti- 
cal lattice. For small and intermediate A, however, a 
magnetization of spin-1 may arise from the triplet states 
on a single lattice site. As a result, the physics of mag- 
netic ordering in double-well lattices at filling two is far 
richer than that in standard optical lattices. 

Magnetic order.- When the (1,1) states dominate the 
on-site Fock states for small A, the interacting Hamilto- 
nian can be mapped to a spin-1 Heisenberg model, which 
can be written as 



H c ff — ^ (JxSr ' Sr+x ~\~ JyS r ■ S r +y) , 



(3) 



where 



2t: 



21- 



2t 2 



2t 2 



u s + u 



Sp 



u v + u, 



sp 



u s + u 



sp 



u v + u s . 



(4) 

S r = A^HA is a spin-1 operator, X is spin-1 Pauli matri- 
ces, and A — (^{, ^J, ^Ij^) 7 are creation operators for 
„t„t in\ JL (J J +p\sl) |0), and 



triplet states p|sj |0), ^= (pjs^ 



Physically, the spin-exchange terms in Eq. ^ come from 
the exchange of fermions with different spins between the 
nearest-neighbor sites in either of the two orbitals. Both 
orbitals contribute to the spin-exchange terms in the ef- 
fective Hamiltonian. 

In the one-dimensional case, Eq. ^ has previously 
been derived in Ref. [TS] For that case, it has been known 
that the one-dimensional spin-1 chain does not have any 
magnetic order, rather the Haldane phase. Nevertheless, 
a two-dimensional spin-1 system can develop antiferro- 
magnetic ground states [2"PH2"2"]. Therefore, one expects 
to see antiferromagnetic spin ordering in a double-well 
optical lattice when the gap A is small and the temper- 
ature is low [23] . 

To characterize the magnetization, we define m = 
— n si + n-pf — n pi\ an d solve the full Hamiltonian 
H = iJband + H[ nt . The results for m as a function of 
the temperature for different A are shown in Fig. |3ja). 
Clearly, the magnetization arises below the Neel temper- 
ature (denoted by T Nco1 ) and saturates to its maximum 
value as the temperature approaches zero. For A = 0, the 
antiferromagnetic ordering is most pronounced: it has 
the highest Neel temperature T Nccl ~ 0.37t. However, 
this must be interpreted with caution because experi- 
mentally it is very difficult to tune the two bands overlap 
with each other while keeping the tight-binding model 
valid. We therefore focus on the cases with nonzero A. 
As A is increased, the magnetization drops faster as T 
increases, and the Neel temperature decreases. For a rel- 
atively large A = 4.5t, T Nccl ~ 0.2U. Note that this 
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FIG. 3: (Color online) (a) Total magnetization m = 
\n s ^ — n s _i + n p f — n P i\ versus temperature, for four selective 
values of A. The chemical potential is selected at approx- 
imate center of the gap (see Fig. |2|. (b) Color plot of the 
magnetization on the A-T plane. In the white regime there is 
no magnetic order, while for low T the magnetization reaches 
the maximum value, indicated as dark colors. The Neel tem- 
perature is shown as blue dashed lines seperating colored and 
white regimes. Note that the y-axis does not start from zero: 
it starts from the lowest temperature T = 0.05t reached in 
the DMFT calculation. 



value of A is already above the critical value A c = At 
in the atomic limit where u(2, 0) — v(0, 2) is the ground 
state. This indicates that many-body effects, such as the 
correlation between nearest-neighbor sites, enhance the 
threshold of A c for a finite m to emerge. 

To give a broader picture, we show in Fig. [3jb) a color 
plot of the magnetization on a plane, of which the axes 
are the energy level splitting A and the temperature 
T. The blue dashed line, separating colored and white 
regimes, indicates the Neel temperature, above which no 
magnetic order is present. The dark color shown near 
the lowest accessible temperature T — 0.05£ in our simu- 
lations characterizes the saturation of the magnetization 
to its maximum value. As the temperature is increased 
to intermediate values, the color turns to red, indicat- 
ing a moderate drop of the magnetization. When the 



temperature is close to the Neel temperature, the mag- 
netization drops rapidly, as can be seen from the narrow 
yellow edge. A close examination of the Neel tempera- 
ture reveals that it drops relatively slowly for A < A c , 
but very rapidly for A > A c . This is consistent with the 
qualitative atomic picture. For A — > oo, the magnetic 
ordering and the corresponding Neel temperature would 
eventually vanish. However, we have shown that the Neel 
ordering would survive at reasonably large values of A. 
This is remarkable, because previous analytical argument 
of mapping to the spin-1 model [IS] is valid for A — » 
only, while for non-zero A a direct numerical solution to 
the lattice model Eqs. Q and Q is highly nontrivial. 
Our results relieve the restriction posed on experiments 
where achieving very small A is difficult. An alterna- 
tive way to appreciate these results is that both m and 
T Nccl are tunable by controlling A, which has obvious 
important experimental implications. 

Enhancement of super- exchange interaction.- The in- 
crease of J [cf. Eqs. gj) and @] in a double- well lattice 
comes from two sources. First, as seen from Eq. Q, in 
addition to t s , t px and t py also enter the expression for 
the super-exchange interaction J. The large value of t px 
then enhances the amplitude of J, similar to the antifcr- 
romagnetism arising from p bands alone [9] . Second and 
most importantly, t sx itself is significantly enhanced in 
a double-well optical lattice. It has been shown that t sx 
can be increased by one order of magnitude at a given lat- 
tice depth for some realistic experimental parameters|14|. 
Thanks to the potential barrier in the center of each lat- 
tice site of a double-well lattice, the Wannier wave func- 
tion of the lowest band spreads its weight toward the 
edge of the corresponding unit cell, which consequently 
enhances the overlap between Wannier wave functions 
on adjacent sites, leading to an increase in the tunnel- 
ing amplitude. As a result, the Neel temperature can 
be strongly enhanced, easily by one order of magnitude. 
The larger energy scale associated with the tunneling and 
super-exchange interaction will also help to reach equilib- 
rim faster in the strongly interacting region. We empha- 
size that this spin-1 antiferromagntism originates from 
the unique feature of the double-well optical lattice: The 
s and p bands can be tuned close to each other, and the 
resulting magnetic ordering incorporates both bands. It 
is this feature that distinguishes our theory from previous 
proposals regarding p bands alone[8l|9]. 

Conclusion.- Using non-perturbative 'DMFT with 
continuous-time quantum impurity solver' direct numer- 
ical techniques, we study two-component fermions in a 
double-well square optical lattice, with two interacting 
orbitals per site. The Mott insulator at filling two is con- 
stituted either by triplet (n s , n p ) — (1, 1) or an admixture 
u(2, 0) — v(0, 2). For the one associated with the triplets, 
antiferromagnetic order emerges in the spin-1 channel be- 
low the Neel temperature, which is determined by the en- 
ergy splitting between the two orbitals and the tunneling 
amplitude. We establish that, as t p contributes to J and 
t sx is significantly enlarged in double-well lattices, the 
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Neel temperature can be one order of magnitude larger 
than that of the one-band system in ordinary optical lat- 
tices, thus perhaps enabling the direct experimental ob- 
servation of the elusive Neel antiferromagnetism in cold 
atomic systems. Our work should facilitate the search of 
magnetic order in optical lattice systems. 
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